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Angle-resolved two-photon photoemission and high-resolution electron energy loss spectroscopy 
are employed to derive the electronic structure of a sub-nanometer atomically precise quasi-one- 
dimensional graphene nanoribbon (GNR) on Au(lll). We resolved occupied and unoccupied elec- 
tronic bands including their dispersion and determined the band gap, which possesses an unexpected 
large value of 5.1 eV. Supported by density functional theory (DFT) calculations for the idealized 
infinite polymer and finite size oligomers an unoccupied non-dispersive electronic state with an en- 
ergetic position in the middle of the band gap of the GNR could be identified. This state resides at 
both ends of the ribbon (end state) and is only found in the finite sized systems, i.e. the oligomers. 

PACS numbers: 81.07.Gf, 73.20.-r, 73.20.Hb, 31.15.-p 



INTRODUCTION 

Along with the great interest in graphene that has 
emerged in recent years other low-dimensional 

carbon-based systems attracted attention, one of which 
are graphene nanoribbons (GNRs). These flat, narrow 
structures are quasi-one-dimensional and have been sub- 
ject to numerous theoretical [2-5] as well as some ex- 
perimental studies [gl-lsl since (unlike graphene itself) 
they exhibit a tunable band gap over a wide range. 
This is important for their possible implementation in 
nanoscale devices such as transistors or transparent elec- 
trodes whereas they maintain the large carrier mobility 
found in graphene sheets [l|,|7|. Besides the width, one of 
the parameters influencing the band gap (and the band 
dispersion along the ribbon axis) is the edge shape. De- 
pending on the corresponding crystallographic direction 
in graphene, the ribbon edge's character varies from arm- 
chair to zigzag quasi-continuously [lOj, much like in car- 
bon nanotubes. The band gap of GNRs are inversely pro- 
portional to the ribbon width 0, [Toj| yielding gap sizes 
up to several electron volts in the sub-nanometer regime. 
Although for a wide range of ribbon widths, this inverse 
relationship has been experimentally observed [f| 0] and 
complies with theory, atomically small widths (<1 nm) 
and the corresponding gaps of more than 1 eV have not 
been observed which is due to the preparation methods 
employed. Another property of interest is the carrier 
mobility in GNRs which manifests itself e.g. in the one- 
dimensional band structure of the ribbons which in first 
approximation resembles the band structure of graphene 
projected onto the respective crystallographic direction 
[2[. However, the lateral confinement of the electrons in 
the ribbon alter the dispersion and particularly open a 
gap at the Dirac point. Theory has dealt in great detail 
with the dispersion in GNRs whereas there have been no 
corresponding experiments published so far. 



Synthesizing suitable GNRs is challenging, particularly 
for the established methods such as lithography [6] and 
unzipping carbon nanotubes [ll| since in both cases the 
widths are usually larger than 10 nm and because of de- 
fects at the edges influencing the electronic structure. 
A powerful alternative in the fabrication of nanostruc- 
tures is the on-surface synthesis [12| , which has been suc- 
cessfully applied for the generation of atomically precise 
GNRs [8]. Thereby the GNR is formed via a surface- 
and thermally assisted [13[ two-step process (see Figure 
□Ja)), in which the adsorbed precursor molecule 10,10'- 
dibromo-9,9'-bianthryl is dehalogenated followed by C-C 
coupling to form a polymer and finally cyclodehydrogena- 
tion yields the defect-free armchair GNRs with a well- 
defined width of 0.7 nm. 

In this letter, we utilize angle-resolved two-photon pho- 
toemission (2PPE) and high resolution electron energy 
loss spectroscopy (HREELS) to determine the electronic 
band structure, i.e. unoccupied and occupied electronic 
states as well as the band gap and dispersion for the GNR 
shown in Figure a) adsorbed on Au(lll). Supported 
by DFT calculations we demonstrate that the band gap 
of the ribbon is surprisingly large namely 5.1 eV and con- 
tradicts all calculated values known from literature [||-[5| . 
An unoccupied electronic state is found in the band gap 
which possesses no dispersion and originates from the 
molecular frontier orbitals localized at the ends of the 
GNR with a finite length. 



RESULTS AND DISCUSSION 

The thermally activated surface-supported formation 
steps of the GNR on Au(lll) can nicely be followed by 
angle-resolved vibrational HREELS, since well-defined 
and pronounced changes in the vibrational spectrum are 
observed for each step. Fig. H{b) shows the data for both 



2 



the polymeric phase and the aromatic nanoribbons. In 




oop edge 
mode 



annealed to T S =470K 
r(c-H) nnear polymer/Au(lll) 

— specular 

— 9°off-specular 



c(C-C) 

I 



;1 



ip-\ 
modes 



v(C-H) 
I 

x 175 



y(c-h) annealed to T S =670K 
GNR/Au(lll) 



oop edge 
mode 



(C-C) 
UJ 



- specular 

- 9°off-specular 



500 1000 1500 2000 2500 3000 

Electron energy loss [cm 1 ] 



FIG. 1. (a) Scheme of the surface-supported fabrication 
of graphene nanoribbons (GNRs) via a two step process: 
(i) dehalogenation of precursor molecule 10,10'-dibromo-9,9'- 
bianthryl followed by a C-C-coupling and (ii) cyclodehydro- 
genation. (b) Changes in the vibrational HREEL spectrum 
observed during the cyclodehydrogenation step from the lin- 
ear polymer to the nanoribbon measured with a primary elec- 
tron energy of 3.5 eV. 

the specular spectrum of the linear polymer, we observe 
several out-of-plane molecular vibrations, namely the C- 
H wagging mode (7(C-H)) at 758 cm -1 , an edge mode 
at 594 cm -1 and the phenyl ring torsion modes (r(C-C)) 
at 424 cm -1 . All these modes show a drastic decrease in 
the spectrum measured in off-specular geometry indicat- 
ing that their intensities are predominantly originating 
from dipole-scattering, i.e., they are dipole active. Thus, 
the corresponding dipole moment changes during vibra- 
tion can be inferred to lie perpendicular to the surface. 
This points towards a mainly flat adsorption geometry 
of the phenyl rings in the polymeric phase. Due to the 
single C-C bonds connecting the building blocks of the 
polymer and the steric hinderance of the C-H groups, 
they are slightly tilted as observed in STM [8]. Indeed, 
at 3055 cm -1 the C-H stretch mode (i/(C— H)) shows a 
weak dipole activity. Accordingly this vibration exhibits 
a component of the dipole moment change parallel to the 
surface normal which is due to a tilting. Hence, our data 
fully supports the adsorption geometry reported by Cai et 
al. Upon cyclodehydrogenation and thus formation 
of the nanoribbons, the vibrational spectrum is changed 
significantly. The number of molecular vibrations found 
in the spectrum is reduced to three modes. The out- 
of-plane modes, the phenyl rings torsion mode r(C-C) 
(404 cm -1 ), the edge mode (609 cm -1 ) and the bend 
mode 7(C-H) at 793 cm -1 are all purely dipole-active. 
This clearly demonstrates that all phenyl rings are now 
orientated parallel (flat-lying) to the surface as expected 



for the aromatic GNR. By adsorption of xenon onto the 
so prepared GNR-covered surface, we may roughly infer 
to the coverage of the nanoribbons on the gold surface: 
in temperature programmed desorption, two desorption 
peaks are observed which do not grow in intensity, as the 
layer thickness is increased to very high coverages. These 
features are assigned to Xe desorbing from both between 
and on top of the nanoribbons, respectively. From their 
peak intensities we can roughly infer to a coverage of 
approximately 2/3 ML. 

Electronic HREELS as well as 2PPE is employed to 
gain insights into the occupied and unoccupied GNR- 
derived electronic states and therefore the band gap. 
Figure [2fa) shows the result obtained from electronic 
HREELS using a primary electron energy of Eo=lbeV. 
A double peak structure is observed which we fitted with 
two Gaussian peaks yielding transition energies of 2.6 eV 
and 3.1 eV, respectively. An exemplary two-color 2PPE 
spectrum recorded with 2.26 eV and 4.50 eV photons is 
displayed in Figure Efb). Several peaks are observed and 
on the basis of photon energy dependent measurements 
[13 . [Til they can be related to photoemission from un- 
occupied intermediate states. The peaks labeled with 
A and B are both populated with the 4.50 eV photons 
and probed by hv=2.26 eV, thus they are located at 4.30 
± 0.06 eV and 3.92 ± 0.06 eV, respectively, with re- 
spect to the Fermi level (Ep)- Note that the energetic 
positions given here are averaged values based on sev- 
eral measurements. The peak labeled C originates from 
a state which is probed by the 4.50 eV photons, there- 
fore it possesses an energetic position of 1.44 ± 0.06 eV 
with respect to Ep. Combining the HREELS transi- 
tion energies with the 2PPE results, we conclude that 
there are two occupied electronic states lying at -1.16 
± 0.08 eV and -1.63 ± 0.08 eV with respect to E F . 
Since these energies do neither comply with the bulk 
band structure of Au(lll) nor with the Shockley surface 
state, these occupied states can be assigned to the GNR. 
Figure [21(c) summarizes the the binding energies of all 
electronic states observed in the present study. On the 
basis of the obtained data one would establish a band 
gap of 2.6 eV as found as the lowest transition in the 
electronic HREELS measurement, which is in reasonable 
agreement with a recent measurement of the band gap 
on spatially aligned GNR using ultraviolet and inverse 
photoemission [16|. However, for an unambiguous identi- 
fication of the band gap one needs further information on 
the nature of the unoccupied bands (states). In order to 
gain deeper insights into the properties we performed dis- 
persion measurements using angle-resolved 2PPE, which 
provide information about the extent of electron delocal- 
ization/localization. 

As observed in STM Q, the GNRs investigated in 
this study are oriented randomly on the Au(lll) sur- 
face. In our experiment, the momentum information is 
obtained by rotating the sample in front of the time-of- 
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FIG. 2. (a) Electronic HREEL spectrum of the GNR-covered 
gold surface recorded at a primary electron energy of 15 eV. 
Two electronic transitions are observed which are fitted using 
two Gaussian peaks, (b) Two-color 2PPE spectra recorded at 
photon energies of 2.25 and 4.50 eV of GNR/Au(lll). (c) En- 
ergetic position of GNR-derived electronic states. The Fermi 
level of Au(lll) serves as reference. 



flight (TOF) spectrometer as seen in Fig. [3fa). While all 
nanoribbons contribute to the signal detected in normal 
emission, this is not the case as we rotate the sample. 
A photoelectron emitted from the one-dimensional band 
structure of GNR band has two contributions to its mo- 
mentum fey parallel to the surface, i.e. the component 
along the ribbon axis (which is the quantity of interest) 
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FIG. 3. (a) Angle-resolved one-color 2PPE spectra (Inset: 
Measuring geometry in angle- resolved 2PPE experiments. 
Among the randomly oriented GNRs on the surface mainly 
those lying in the plane perpendicular to the rotation axis of 
the sample are detected for off-normal emission), (b): Dis- 
persion of the spectral feature labeled as A, B and C. While 
the states A and C are localized, state B shows a strong dis- 
persion, viz. is delocalized. 



and a component perpendicular to it which is random (as 
the electrons are localized in this direction). Thanks to 
the measuring geometry our TOF only detects photoelec- 
trons emitted in the plane perpendicular to the rotation 
axis. This can either be electrons emitted along the axis 
of ribbons lying in this plane or electrons which (due to 
the random momentum perpendicular to the ribbon axis) 
are emitted at an angle with respect to the ribbon which 
is equal to the angle between this GNR and the detection 
plane. In the latter case, the detected momentum will be 
higher than the momentum the electron had in the band 
which causes the measured dispersion to be smeared out 
to higher k but not in energy. Since for finite angles only 
ribbons lying perpendicular to the rotation axis can be 
measured, we expect a significant drop in photoemission 
intensity away from the F point. 

From the three unoccupied bands observed in 2PPE 
only the one labeled as B shows a dispersion, viz. the 
states A and C are localized. Figure [3fa) displays a series 
of one-color 2PPE spectra recorded at a photon energy 
of 4.66 eV for various angles between the surface nor- 
mal and the TOF axis. Beside the contributions of the 
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gold d-bands, the peaks labeled A, B, and C are seen. 
While B exhibits a strong dispersion around the T point, 
state A and C show no dispersion (see Fig. [3fb)). The 
parabolic behavior of state B complies qualitatively with 
all theoretical calculations on armchair nanoribbons. The 
effective mass (m*) of 1.2 m e is significantly higher than 
theoretical models generally predict for a variety of GNRs 



state C could be found in the periodic calculations for 
the idealized infinite polymer. Therefore, we carried 
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However most surprising the lowest unoccupied state 
located at 1.44 eV above Ep (peak C) exhibits no disper- 
sion, therefore we assign the first dispersive state (state 
B) above the Fermi level, which possesses an energetic 
position of 3.92 eV to the conduction band of the GNR. 
Thus, with the valence band located at -1.16 eV with 
respect to E F the band gap is 5.1 eV. This value is sig- 
nificantly higher than the calculated values lying in the 
range between 1.5 eV and 2.1 eV depending on the ap- 
plied method [8]. For the use in technological applica- 
tions, the size of the band gap plays a key role; the one 
we observe for the studied armchair GNR is too large for 
instance the use in a nanowire or even a transistor where 
a small gap is preferred. 

In order to understand the nature of the localized state 
C we performed DFT calculations for the idealized infi- 
nite polymer and finite size oligomers employing the PBE 
GGA-functional and the range separated hybrid func- 
tional CAM-B3LYP, which is particularly suitable for 
long chain molecular systems (polyene) [19] together with 
the 6-31 1G** basis set as implemented in the Gaussian 09 
program package [20[ . All calculations were done for the 
GNR alone, i.e. without the supporting Au(lll) surface, 
and the geometry of the GNR was fully optimized. For 
the periodic calculations a /c-point grid with 240 points 
in the first Brillouin zone (1BZ) and a real space cut off 
of 600 bohrs were used. The PBE results obtained here 
compare very well with calculations reported by Cai and 
coworkers [8[- Furthermore, for the periodic calculations, 
the geometry of the GNR depends only very weakly on 
the choice of the functional. For instance, we obtain a lat- 
tice constant of 4.310 A for PBE and for CAM-B3LYP we 
get 4.277 A. However, the band gap changes drastically 
due to the admixture of exact exchange in the CAM- 
B3LYP functional. Using the PBE functional results in 
a gap of 1.54 eV, while for CAM-B3LYP a gap of 4.35 eV 
is obtained. The same qualitative trend is observed when 
using PBEO [2l[, i- e -? PBE with an admixture of exact 
exchange, where we obtain a band gab of 2.7 eV. Also, it 
has been shown, that CAM-B3LYP is especially suitable 
for structural and electronic properties of large unsatu- 
rated hydrocarbons [22[, which is in line with the much 
better agreement of CAM-B3LYP with the experimental 
results reported here. 

However, apart from the wider band gap and the 
larger dispersion, the band structures for the PBE and 
CAM-B3LYP are qualitatively similar and no indication 
of a non-dispersive state corresponding to the localized 
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FIG. 4. (a) Density of state obtained be Lorentzian- 
broadening for the periodic case (black line) and for the oc- 
cupied (blue lines) and virtual orbitals (red lines) of the hep- 
tamer and the hexamer together with the experimental ener- 
gies, which have been aligned to the middle of the band gap 
of the periodic GNR. (b) HOMO and LUMO for the hep- 
tamer, which are localized at the ends of the oligomer (all 
calculations at the CAM-B3LYP/6-311G** level of theory). 

out additional calculations for oligomers up to the hep- 
tamer, where the oligomers were saturated with hydro- 
gens at both ends and again a full geometry optimiza- 
tion was performed. In order to compare the electronic 
structure of the oligomers with the one of the periodic 
structure, we computed the density of states (DOS) by 
a simple Lorentzian-broadening procedure, where every 
Kohn-Sham orbital energy or every crystal orbital en- 
ergy within the 1BZ is represented by a Lorentzian with 
a width T=0.2 eV for the oligomers and T=0.1 eV for 
the periodic case. Figure IHa) shows the DOS for the pe- 
riodic case (black line) and for the occupied (blue lines) 
and virtual, i.e., unoccupied, orbitals (red lines) of the 
heptamer and the hexamer using the CAM-B3LYP func- 
tional. As one can see, the DOS for the valence and the 
conduction bands are quite similar. The same holds for 
the geometric structure, as we get a length of 4.278 A for 
the monomer unit in the middle of the heptamer, which 
is nearly identical to the lattice constant of 4.277 A in 
the periodic case. 

However, for the electronic structure of the oligomers 
we observe two states inside the "band gap" , namely the 
highest occupied (HOMO) and lowest unoccupied molec- 
ular orbital (LUMO). These states are localized at both 
ends of the oligomers, as shown in Fig. SJb) for the hep- 
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tamer. Furthermore, the energy difference between these 
two localized states decreases for longer chain lengths. 
In the experiments GNRs with a length, which corre- 
sponds to dozens of units (monomer length: 0.43 nm) 
[8[, for which a non-resolvable energy difference between 
the two localized states is to be expected. In addition, 
the HOMO of the free GNR should also be unoccupied, 
when adsorbed at the Au(lll) surface, as its energy is 
above the Fermi level. Therefore, these localized states 
can be identified as state C comparing the calculated en- 
ergies with experimental determined energies, which are 
indicated at the right side of Fig. [4f a). 

CONCLUSION 

In conclusion, using two complementary surface- 
sensitive experimental methods, namely HREELS and 
2PPE, in combination with DFT calculations we were 
able to probe key electronic properties of a defect-free 
graphene nanoribbon. We determined the band gap to 
be 5.1 eV, which is surprisingly high since previous cal- 
culations predicted much lower values. A non-dispersive 
unoccupied electronic state of the GNR located in the 
band gap originates from both the HOMO and LUMO 
of the ribbons with a finite length (oligomer). These 
states are localized at both ends of the ribbon, thus we 
named them "end states". A controlled modification of 
the precursor molecule with suitable substituents or by 
doping may certainly open the perspective to create pre- 
cise nanoribbons with lower band gaps, appropriate for 
applications. 
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